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The d + N systems are studied in a three-body model, using phenomenological N-N interactions. 
The scattering matrices are calculated by using the Kohn-Hulthen variational method. Then, they 
are analytically continued to complex energies and their singularities are localized. We find a virtual 
state at E = —1.66 MeV in and a pair of states at E = (—0.42 ± *0.52) MeV in ®He relative to 
the d + N thresholds, respectively. In addition, we discuss some general aspects and problems of 
virtual states which may be useful also in the study of other systems, like ^°Li. 
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I. INTRODUCTION 

The lightest nucleus with a well-established spectrum 
of excited states is '*He. There have been numerous at¬ 
tempts, both experimental and theoretical, to identify 
excited states in A = 3 nuclei; for an extensive review see 
Ref. [|p . Although there are indications for resonance-like 
behavior in the 3n and 3p systems in certain experiments 
1^ , others do not seem to support these findings . The¬ 
oretically the most well-established case is a broad three- 
neutron (and three-proton) resonance in the J'^ = 3/2+ 
partial wave jij. However, the existence of these struc¬ 
tures has not been confirmed by experiments yet. 

According to Ref. Q, no unambiguous experimental 
evidence has been found to support the existence of low- 
lying and ^He resonances. However, several calcula¬ 
tions, mostly conhned to small model spaces, indicate a 
virtual state in ^H, and corresponding subthreshold res¬ 
onances in ^He. Recently, the authors of a H(®He, "^He) 
experiment have claimed to find an excited state of ^H at 
A* = 7 ± 0.3 MeV excitation energy with T = 0.6 ± 0.3 
MeV width . According to one interpretation this is 
a 1/2+ excited state formed mainly in the (rin)p channel. 
In the mirror ^He nucleus the Ref. experiment does 
not indicate an excited state. However, in a recent work 
0 it was shown that the effective range function of the 
doublet d-\-p scattering has a singularity at a negative en¬ 
ergy, just below the d-\-p threshold. This would indicate 
the existence of subthreshold virtual states in ^He. 

In the present work we study the ^H and ^He nuclei 
with the aim to clarify some of the possibilities men¬ 
tioned above. We describe the d N scattering in a 
nonrelativistic three-body potential model and study the 
analytic properties of the scattering matrices. Since vir¬ 
tual states play an important role in our study, we briefly 
discuss some of their features in Sec. H, and show some 
problems associated with their experimental investiga¬ 
tion. Our three-body calculations are reviewed in Sec. 
HI, while Sec. IV offers some conclusions. 


II. VIRTUAL STATES IN SCATTERING THEORY 

The spectrum of a quantum scattering system gener¬ 
ally consists of bound states and continuum states. The 
continuum has structure caused by the complex-energy 
singularities of the scattering matrix. For a broad class of 
physically interesting interactions all singularities within 
the meromorphic region of the potential are physical and 
correspond to solutions of the Schrodinger equation with 
purely outgoing asymptotics. The following types of sin¬ 
gularities are usually defined on the complex k (wave 
number) plane: i) bound states with ks = iy (7 > 0); 
ii) virtual (antibound) states with kA = —*7 (7 > 0); 
resonances with kn = — 17 (k, 7 > 0; k > 7); and 

quasi-resonances with kc = ±«: — iy (k, 7 > 0; k < 7). 

Our main interest in the present work is in the vir¬ 
tual states. Interestingly, although a virtual state is an 
unbound state with exponentially increasing wave func¬ 
tion, it corresponds to negative energy, Ey = k’^h?/2m = 
—7^?i^/2m. The cross section of a process which involves 
a virtual state behaves like cr{E) ~ 1/(A -|- \Ev\) for 
E > 0. Thus it is singular at the unphysical negative en¬ 
ergy E = —\Ev\, and increases with decreasing positive 
energies. This behavior can have rather strange effects 
on experimental measurements. If a measurement is per¬ 
formed on a system which has a virtual state and the 
interaction-free cross section (phase space) of the mea¬ 
sured process drops to zero with decreasing energy, then 
the result of the measurement is a cross section with 
a low-energy peak. This may be interpreted as a low- 
energy resonance at positive energies. However, the peak 
cannot be fitted with a Breit-Wigner form, which causes 
the extracted “resonance parameters” to be highly de¬ 
pendent on the fit procedure. A good example of this 
experimental difficulty concerning virtual states is ^°Li. 
We believe that the experimental confusion regarding the 
nature of the ground state of ^°Li B is caused by the fact 
that there is a virtual state in the ^i-|-n system |^. The 
energy of this state is negative, for instance Ey « —0.03 
MeV for the P2 interaction of Ref. Q. If, however, the 
experimental procedures assume a low-energy resonance 
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in this system and try to fit the various measured cross 
sections accordingly by Breit-Wigner forms, then they 
can get rather different results. 

We would like to emphasize that pure virtual states 
with pure imaginary wave number can only exist in neu¬ 
tral systems with two-body asymptotics and only in S- 
wave states. By making the interaction more attrac¬ 
tive the virtual state can be moved from the negative 
imaginary fc-axis to the positive one, ending up with a 
bound state [0. If, however, a Coulomb or centrifu¬ 
gal barrier is added to the interaction, then the virtual 
pole moves from the imaginary axis into the complex 
k plane resulting in a quasi-resonance. In order to il¬ 
lustrate this process we show in Fig. the pole tra¬ 
jectory corresponding to the ^Sq N + N system. We 
start from n + n and smoothly switch on the p + p 
Coulomb interaction to finally end up as the p + p sys¬ 
tem. One can see that the effect of the Coulomb force 
is to move the virtual pole into the complex plane (and 
to create a conjugate pole) as discussed above. We use 
here the Eikemeier-Hackenbroich (EH) N-N interaction 
[pd| . The n + n and p + p poles are at Enn = —0.134 
MeV and Epp = (—0.101 ± i0.515) MeV energies re¬ 
spectively. These numbers are in reasonable agreement 
with those that can be extracted from experimental data, 
= -0.123 MeV and Epp = (-0.140 ± fO.467) MeV 
[p^ . Our interaction is charge-independent, so it gives 
the same pole position for n + n and n + p. Naturally, it 
cannot reproduce the experimental Enp = —0.066 MeV 
0 - 

We note that although one cannot have pure virtual 
poles in charged systems, a pair of quasi-resonance poles 
close to the negative imaginary k axis can have observable 
effects at positive real energies, like in the case oi p + p. 
The reason for this is the fact that unlike in the case of 
usual resonance poles, both of such poles are roughly the 
same distance from the real energy axis of the physical 
energy sheet. This is not the case for usual resonant 
states. The e = Er — iT/2 {Er,T > 0) main pole of a 
usual resonance can be reached by simply crossing the 
real energy axis and going to the fourth quadrant of the 
energy plane. However, the conjugate pole can only be 
reached through a long path by crossing the real energy 
axis, going into the fourth quadrant, and then encircling 
the origin and ending up at e = -I- ir/2 {Er, T > 0). 

To recap our brief discussion of virtual states, we em¬ 
phasize the sometimes forgotten fact that the energy cor¬ 
responding to these states is negative. We also note that 
in charged systems pure virtual states cannot exist. How¬ 
ever, the effect of a pair of conjugate quasi-resonances 
close to the negative imaginary k axis (corresponding to 
a virtual state in a neutral system) can have significant 
observable effects. 


III. THREE-BODY CALCULATIONS FOR d-|-N 
SCATTERING 


We solve the three-body Schrodinger equation for ^H 
and ^He using the EH N-N interaction. This interaction 
gives a good general description of the N + N scattering 
data and the deuteron properties |l|]. Technical prob¬ 
lems prevent us from using the most modern interactions, 
but for our current purposes the EH interaction is quite 
satisfactory. In our calculations the same nucleon mass 
Mjv = {Mp + M„)/2 is used for both protons and neu¬ 
trons. 

In this paper we study only = 1/2+ states. As a 
first step we solve the bound-state problem of ^H and ^He 
using the variational Gaussian-basis coupled rearrange¬ 
ment channel method HI . We include all 23 channels 


with li,l 2 < 2 in the 


(^ 1 , 52 )^ 12 , ^3 S,{h,l2)L 


coupling scheme. Here Si = S 2 = S 3 = 1/2 are the 
nucleon spins, S 12 is the coupled two-nucleon intrinsic 
spin, S is the the total intrinsic spin, li and I 2 are the 
orbital angular momenta of the two relative motions, 
respectively, L is the total orbital angular momentum, 
J is the total spin, and tt = (—l)*i +*2 jg parity. 
Our calculated binding energies, Es^ = —7.65 MeV and 
i?3He = —6.99 MeV, are smaller than the experimental 
values = -8.482 MeV and = -7.718 MeV, 
respectively. Our triton binding energy is close to those 
coming from calculations that use the most modern two- 
nucleon forces Thus, a large part of the underbind¬ 
ing can be attributed to the lack of three-body forces in 
our model. 

We find that by keeping only the 9 most important 
channels, shown in Table |, the binding energies change 
by only 0.05 MeV. In the following we keep only these 9 
channels. In order to have a coupled — ^Di deuteron 
in our asymptotic d + iV-type wave function compo¬ 


nents, we use the 


(^1, ^2)^12, 


h,S3 

I,l2 




cou¬ 


pling scheme in the scattering calculations. Here Ii is 
the total (intrinsic plus orbital) spin of the two-nucleon 
subsystem and I comes from the coupling of /i and S 3 . 
We show our model space in this coupling in Table ||. 
Hence, we have 7 channels: two channels containing a 
— ^Di deuteron plus a nucleon with I 2 = 0 and 2, 
respectively (lines 1-2 and 3-4 in Table ||), three chan¬ 
nels with states inside the n + p system (lines 5,6, 
and 7 in Table I® , and two channels with ^ Sq n+p and 
n + n {p + p for^He) two-body subsystems, respectively 
(lines 8 and 9 in Table ||). The multichannel scattering 
problem is solved by using the Kohn-Hulthen variational 
method for the S matrix El- In all our calculations 
we remain below the three-body breakup threshold. It 
means that in channels 5-9 in Table || the wave func¬ 
tions have three-body bound-state asymptotics |^. In 
these channels we keep the same Gaussian basis as in the 
bound state calculations. 
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The scattering matrices, coming from the Kohn- 
Hulthen calculations, are analytically continued to the 
multisheeted complex Riemann energy surface using the 
methods of Ref. [|^, and their singularities are local¬ 
ized. As a test calculation we searched for the bound- 
state poles and found them at the same energies as in 
the bound-state calculations. In addition, a virtual state 
is found in at = —1.66 MeV, relative to the d + n 
threshold. We note that as our triton ground state is 
underbound so is probably this virtual state. The ef¬ 
fect of underbinding in a virtual state is that \Ey\ be¬ 
comes larger 1^. This means that if our N — N in¬ 
teraction reproduced the correct binding energy of 
and ^He then the virtual-state pole would be closer to 
zero energy. In ^He we find a pair of conjugate poles at 
= (—0.42±i0.52) MeV, relative to the d-\-p thresh¬ 
old. The insufficient attraction in our model probably 
moves the poles too far away from the negative energy 
axis. 

A practical way to extract the parameters of the vir¬ 
tual states of and ^He from experimental data is 
to describe these nuclei within the extended i?-matrix 
model [l^ ], which can be continued to complex energies. 
An extensive analysis of the A = 3 data is being per¬ 
formed using this method M]. Preliminary results indi¬ 
cate the presence of and ^He virtual states at 

= -1.07 MeV and A;^^® = (-0.72 ± iO.23) MeV, 
respectively. The results of our three-body calculations 
are consistent with these findings if we take into account 
the fact that the insufficient binding probably increases 
\E^\ and pushes the ^He pole away from the negative 
axis. 

In comparing our results to previous calculations, we 
note that in Ref. [^, negative-energy poles were found 
in the ^S' 1/2 kcotd function for d -I- n at w —160 keV 
and for d -I- p at « —25 keV, using the S-wave MT I- 
III interaction of Ref. Ill- A similar pole was found in 
the d -I- p phase-shift analysis of Black et al. . These 
real-energy poles in the d + N effective-range expansion 
are clearly related to the presence of the virtual states 
we have found, as will be discussed in a later publication 

0 . 

IV. CONCLUSION 

In summary, we have localized virtual d -I- V states in 
and ^He in a model that contains all relevant three- 
body channels. These states can be considered as special 
excited states of the A = 3 nuclei. Our results undoubt¬ 
edly show the existence of these structures, although 
their properties could change somewhat if the most mod¬ 
ern N-N interactions with three-body forces were used. 
Such calculations are, however, beyond the scope of the 
present work. 
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TABLE 1. Channel configurations used in our 9-channel 
calculations in the [(Si, S'2)5'i2, Ssj S', (L,/ 2 )L coupling 


scheme. 


Configuration 

Sl2 

S 

h 

h 

L 

N{pn) 

1 

1/2 

0 

0 

0 

N{pn) 

1 

3/2 

2 

0 

2 

N{pn) 

1 

3/2 

0 

2 

2 

N(pn) 

1 

1/2 

2 

2 

0 

N{pn) 

1 

1/2 

2 

2 

1 

N(pn) 

1 

3/2 

2 

2 

1 

N(pn) 

1 

3/2 

2 

2 

2 

N{pn) 

0 

1/2 

0 

0 

0 

p{nn) or n(pp) 

0 

1/2 

0 

0 

0 


TABLE 11. Channel configurations 

nel calculations in the (Si, S 2 )Si 2 , L 

pling scheme. 

used in 

. /l,S3 

our 

/,/2 

9-chan- 

cou- 

Configuration 

Sl2 

h 

h 

I 

h 

N{pn) 

1 

0 

1 

1/2 

0 

N{pn) 

1 

2 

1 

1/2 

0 

N{pn) 

1 

0 

1 

3/2 

2 

N{pn) 

1 

2 

1 

3/2 

2 

N{pn) 

1 

2 

2 

3/2 

2 

N(pn) 

1 

2 

2 

5/2 

2 

N(pn) 

1 

2 

3 

5/2 

2 

N{pn) 

0 

0 

0 

1/2 

0 

p{nn) or n{pp) 

0 

0 

0 

1/2 

0 



- 0.1 0.0 0.1 


Re(k) (fni') 

FIG. 1. Trajectories of the ^Sq N + N S-matrix poles. The 
open square corresponds to the n-\-n and n+p poles, wile the 
open circles denote the pair of conjugate poles in the p-\-p sys¬ 
tem. The filled circles come from calculations where c • 
is added to the n -|- n interaction (0 < c < 1). 
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